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ABSTRACT 


This report discusses a new and innovative numerical procedure for evaluating and 
monitoring contact stresses in meshing gear teeth. The procedure is intended to extend the 
range of applicability and to improve the accuracy, of gear contact stress analysis. 

Originally proposed by Burton Paul, et al., the procedure is based upon fundamental 
solution from the theory of elasticity. It is an iterative numerical procedure. The method 
is believed to have distinct advantages over the classical Hertz method, the finite-element 
method (FEM), and over existing approaches with the boundary element method (BEM). 

Unlike many classical contact stress analyses, friction effects and sliding are included. 
Slipping and sticking in the contact region are studied. 

Several examples are discussed. The results are in agreement with classical results. 

Applications are presented for spur gears. 
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1. INTRODUCTION 


Tooth contact stresses are not nearly as well understood as root stresses and fillet 
stresses, and yet contact stresses are among the most important factors in determining gear 
life. Contact stresses also have a significant effect upon tooth deformations and thus upon 
gear performance under load. 

One of the reasons that gear contact stresses are not better understood is the 
geometry and kinetics of meshing gear teeth do not satisfy the assumption required of the 
direct application of established Hertzian procedures. The Hertz analysis assumes contact 
between frictionless cylinders. Also, since contact stresses represent stress concentrations 
with nonlinear behavior, accurate results are not readily obtained with standard finite 
element procedures. Finally, the use of other numerical procedures such as the boundary 
element method (BEM) are not yet fully developed so that they are easily applied with 
mating gear tooth surfaces. 

In this report we present a new contact stress analysis procedure which is intended 
to produce an improved understanding of gear contact stresses. This procedure does not 
require the restrictive assumptions of the Hertz analysis. Specifically, with this new 
procedure non-cylindrical geometry can readily be accommodated and frictional effects can 
be included. 

The approach is based upon the point load superposition method developed by Paul, 
et. al. [21,24,27] and previous, unpublished work of Sisira Jayasinghe at the University of 
Cincinnati. The report itself is an expansion of a paper presented at the NASA sponsored 
Second Annual Conference on Health Monitoring for Propulsion Systems in Cincinnati 
(1990) [47], 

Historically, contact stresses have been one of the most important, most widely 
studied, and yet most difficult problems in the theory of elasticity. Contact stress analysis 
can be traced back to the work of H. Hertz in 1881 [1,2,3] who studied the contact 
deformation of smooth elastic bodies pressed together by forces normal to their contacting 
surfaces. He proposed that the contact region could be modelled by an elastic half space 
loaded over an elliptical area. 

Hertz’s theories have formed the basis for contact stress analyses for nearly a century. 
However, important contributions, extensions, and applications have been developed by 
Ahmadi, Keer, and Mura [4]; Barber [5]; Bentall and Johnson [6]; Boussinesq [7]; Bryant 
and Keer [8]; Carter [9]; Cerruti [10]; Conry and Seireg [11]; Cooper [12]; Danders [13]; 
Dyson [14]; Gladwell [15]; Goodman [16]; Johnson [17, 18]; Kalker and Piotrowski [19]; 
Keer, Dunders, and Tsai [20]; Liu and Paul [21]; Ix> [22]; Mindlin [23]; Paul and Hashuni 
[24]; Paritsky [25]; Sackfield and Hills [26]; and Singh and Paul [27]. Excellent summaries 
of contact stress analysis can be found in the books of Gladwell [28] and Johnson [29] and 
in the review article of Kalker [30], 

References [31 to 44] describe recent attempts to use the finite element method 
(FEM) and the boundary element method (BEM) to determine contact stresses. Since 


’Numbers in brackets refer to the References. 
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contact stresses form local stress concentrations, there are convergence problems with the 
FEM due to the large stress gradients. Many elements are needed in the contact region to 
obtain an accurate representation of the stress. Also, iterative procedure and nonlinear gap 
elements are needed to define the contact region. Alternatively, the BEM, which is useful 
for problems with high stress gradients, also requires an iterative procedure to determine 
the contact region, but it may provide more accurate results than the FEM. Unfortunately, 
there is relatively little BEM software, with documentation, applicable to contact analysis. 

The method of this report is an iterative numerical procedure based upon 
fundamental solutions of the theory of elasticity. As with the BEM it employs boundary 
elements on the surfaces of the mating bodies, but only within the contact region itself. 

The balance of the report is divided into nine parts with the following part 
summarizing some preliminary results useful in the sequel. The next three parts present the 
basic analysis. Numerical procedures are then discussed followed by some simple examples. 
The next part discusses application with spur gears. The final two parts contain a discussion 
and concluding remarks. 


2. PRELIMINARY CONSIDERATIONS 

Consider two elastic bodies in contact and let the contact region be divided into a mesh 
of elements (or cells). Let the forces transmitted across each cell be represented by a 
distribution of concentrated forces acting on an elastic half space. Then by using classical 
results from the theory of elasticity, together with the principle of superposition, the stress 
distribution and point displacements due to a cell loading are obtained. Finally, by 
superposing the results from each cell the resulting stress and displacement distributions 
within the contacting bodies are obtained. 

To develop the analysis consider an elastic half space with a concentrated normal load 
as in Figure 1. 



Figure 1. A Concentrated Normal Force on an Elastic Half Space 
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Let P be a typical point within the half space, located by the polar coordinates (r, 0) 
as shown. Then the radial and tangential stresses at P are [29, 45, 46]: 

2F 

o rr cos 0, o re - o e0 - 0 (1) 

TZ T 

In rectangular coordinates, these stresses have the form: 

2F x 2 z 2F z 3 2F xz 2 

a - - — , o - - — , o (2) 

“ * (x 2 + z 2 ) 2 “ * (x 2 +y 2 ) 2 “ * (x 2 +z 2 ) 2 

Using the stress-strain equations and the strain-displacement equations, the 
displacements at a point P on the surface may be expressed as [29]: 

_ p (l-2v)(l+v) . 2F to!h/xl (3) 

* 2E |x| icE 

where E and v are the elastic constant and Poisson’s ratio, and where h is the distance to 
a reference point of zero displacement along the z-axis. 

In like manner, if there is a concentrated tangential force G acting at o the stresses at 
P are [29]: 

o rr - - — cos 0, o fl - o - 0 (4) 

rcr 


where 0 is now measured from the line of action of G. In rectangular coordinates these 
stresses are 


2G x 3 2G xz 2 2G x 2 z 

* (x 2 +z 2 ) 2 ’ ' " * (x 2 +z 2 ) 2 ' " * (x 2 +z 2 ) 2 


(5) 


The displacement of a surface point are then 


(1-v 2 ) . | | n (l-2v)(l+v) . x . 

u. - -2G to lx I + c, u - G — (-r-r) 

x nE 2 2E Ixl 


(6) 


where c is a constant determined from a zero-displacement reference point 

Next, consider an elastic half space with a distributed normal and tangential loading over 
a portion of the surface as depicted in Figure 2. 
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Figure 2. Elastic Half-Space with a Normal and Tangential Loading Distribution 

Then the resulting stresses and displacements at a typical point P may be obtained by 
superposition (integration) of the foregoing results. Following the analysis of Johnson [29], 
the stresses at P are: 

- - — f[ {[zf(s)(x-s) 2 + g(s)(x-s) 3 ]/[(x-s) 2 + z 2 ] 2 )ds 

It J ~ b 


a zz - - — [[ l[z 3 f(s) + z 2 g(s)(x-s)]/[(x-s) 2 + z 2 ] 2 lds 

It J b 

a xz — f 1 {[z 2 f(s)(x-s) + zg(s)(x-s) 2 ]/[(x-s) 2 + z 2 ] 2 )ds (7) 

n J * b 


In the manner the tangential and normal displacements of a surface point are [29]: 

u - - (1 ~ 2v)(1+v) [f x f(s)ds - f* f(s)ds] - 2(1 ~ v2) - f* g(s)to|x-s|ds+c, 

x 2E J- b Jx nE J ~ b 

( 8 ) 

. . /• f( s)to |,- s |ds . g(s)ds - /• g(s)ds]+c 2 

7tE 2E Jb Jx 
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where Cj and c 2 are constants to be determined from the displacement of an arbitrary 
reference point. 

3. APPLICATION: TRIANGULAR SHAPE LOADING DISTRIBUTION 

Equations (7) and (8) may be used to determine the effects of the force distribution in 
the contact region. If the region is divided into a mesh of elements (or cells), then the 
effects of the loadings on the individual cells may be superimposed. Moreover, the general 
loading in contact region may be represented by the superimposed loadings on the cells. 
A triangular (that is, piecewise linear) loading distribution over a cell forms a convenient 
basis for modelling the general loading distribution. 

The contribution of a triangular load on an individual cell to the stress at a typical point 
P may be obtained by letting f(x) and g(x) of Equations (7) and (8) have the forms: 

f(x) - f 0 (a-lxl)/a , g(x) - g,j(a- lx l)/a , Ixl s a (9) 


where a is the half-width of the cell whose center is at the origin, and where f 0 and go are 
the peak triangular normal and tangential loadings as depicted in Figure 3. 



Figure 3. Triangular Normal and Tangential Loading on a Topical Cell 

By substituting from Equations (9) into (7) and by performing the indicated integrations, 
the resulting stresses at P (see Figure 3) are [29]: 
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( 10 ) 


f 

a ^ [(x-a)0 1 + (x+a)0 2 - 2x0 + 2zln(rjr 2 /r 2 )] 

Tta 

+ — [(2xln(r 1 r 2 /r 2 ) + 2aln(r 2 /rj) - 3z(0, + 0 2 - 20)] 

Tta 

- — [(x-a)0 1 + (x+a)0 2 - 2x0] - -^(0, + 0 2 - 20) 
ixa ita 

o “ — ——(0. + © 2 — 20) + — [(x-a)0j + (x+a)0 2 - 2x0 + 2zln(rjr 2 r 2 )] 
rca rca 


where r 1; r 2 , r, 0^02, and 0 are given by (See Figure 3): 
rf - (x-a) 2 , r 2 - (x+a) 2 + z 2 , r 2 - x 2 + z 2 

( 11 ) 

tan0j - z/(x-a) , tan0 2 - z/(x+a) , tan0 - z/x 

Similarly, the displacements of a surface point are: 
u K - Ag 0 [(x+a) 2 fo(x/a+ 1) 2 + (x-a) 2 to(x/a-l) 2 - 2x 2 to(x/a) 2 ] + Bf 0 + C t 

( 12 ) 

u z - Af 0 [(x+a) 2 to(x/a+l) 2 + (x-a) 2 Cn(x/a-l) 2 - 2x 2 to(x/a) 2 ] - BgQ + C 2 

where Q and Q are constants to be determined from the displacements of a reference point 
and where A and B are defined as: 
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a -dk*i 

2nEA 


B 


' ^ (l-2v)(l+v) a 
< E 2 

_ (l-2v)(l+v) x ( a - 1 a/2 | ) 
Ea 


for x ^ a 
for | x | ^ a 


(13) 


4. MODELLING OF CONTACT FORCE DISTRIBUTION 

To illustrate how these results can be superimposed in the contact region, suppose that 
the contact force distribution is represented by the piecewise linear distribution, as in 
Figure 4. 



Figure 4. A Piecewise Linear Representation of a Force Distribution 

Let the contact region be divided into n equal-width elements (or cells) as shown. 
Then, just as the region is discretized so also are the force and displacement distributions. 
Let a be the element half-width. Then the element nodal coordinates are in multiples of 
a. That is, the x-coordinate of node k is: 

x k - ka (14) 

Let u^ and u^ represent the tangential and normal displacements of node i due to 
a triangular normal load centered at node j. Then from Equations (12), u jx and u, z may be 
expressed as 


t 
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(15) 


_ _ (l-2v)(l+y) d f 

IX £ IJ J ’ 


u ; . 


.l- v2 v f 
- <— 


— [(k+ l) 2 ln(k+ l) 2 +(k- l) 2 ln(k- l) 2 -2k 2 lnk 2 ] 
2n 

&d>l 


0 

k-0 

(16) 

2a .. 

to4 

. n 

k-±l 



' a/2, k>0 
< -a/2, k<0 ► 
0, k-0 


where k is defined as: 
k - i - j 


(17) 


Equation (15) may be used to obtain the displacements at node i due to a set of 
superposed triangular load distributions centered at each of the nodes. In this case the 
terms c,^ and d if fj are interpreted as sums over j. 

Finally, if there is a horizontal or tangential loading, an analogous analysis leads to 
displacements at node i as: 


u ix - (-^-) C ijgj 


u. 


(l-2v)(l+v) 

E ij j 


(18) 


where the coefficients c^ and dy are the same as those of Equations (16) and where there 
is a sum on j from 1 to n. 

The displacements of Equations (15) and (18) are measured relating to a convenient 
reference point, say the origin of the coordinate axes. 

5. CONTACT ANALYSIS 

Equation (15) and (18) may be used to develop an iterative procedure for obtaining a 
detailed and comprehensive analysis of contact stresses and deformations. 
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Figure 5. Two Elastic Cylinders Being Pressed Together 

To illustrate this consider two cylinders pressed together as in Figure 5. As the cylinders 
are brought together they initially have line contact. Then as they are pressed tighter the 
line contact develops into a contact strip. 

Figure 6 depicts the contact region after a finite deformation. The solid curves 
represent the deformed surface profiles, and the dashed curves represent the undeformed, 
or original profiles. 



Figure 6. Contact Region of Deformed Cylinders 
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Consider two points Qj and Q 2 relatively distant from the contact region. Then during 
the compression Qj and Q 2 move toward the contact region and toward each other through 
the displacements d lz and d 2z respectively. 

Consider two typical matching surface points P x and P 2 separated by a distance h(x). 
Then as the deformation proceeds, P, and P 2 simultaneously move toward each other and 
deform inward into their respective cylinders. Let U lz and U 2z (measured as positive into 
the respective cylinders) represent the z-axis displacements of Pj and P 2 toward Q, and Q 2 . 
Then if P t and P 2 come into contact, we have 

U„ * U 2l + h(x) - d„ * d* - D, (19) 

where D z is the global compressive displacement of the cylinders. If P t and P 2 are outside 
the contact region we have 

U lz + U 2z + h(x) > D z (20) 

Equations (19) and (20) can be used to determine the extent of the contact region. 

If there is no relative horizontal displacement of P t and P 2 , we have the no-slip 
condition 

U u - u 2l < 21 > 

When the cylinders are subjected to tangential loading in addition to the normal loading, 
the deformation geometry depends upon the relative magnitudes of the normal and 
tangential loads. The tangential loading causes on asymmetric deformation pattern. This 
asymmetry, or distortion is limited by the friction coefficient so that alternatively sliding 
occurs. 

Figure 7 depicts the cylinders in contact with tangential loading. As before ?! and P 2 
are matching surface points at the initiation of contact. When the tangential load G is 
applied the distant points Qj and 0 2 have "rigid body" displacements d lx and d 2x relative to 
O. If the displacements of Pj and P 2 relative to O are denoted by S lx and S 2x , the "slip" S x 
between P t and P 2 may be defined as: 

S„ - S„ - s 2i - (U„ - d, ( ) - (U 2 , - <y - (U„ - u 2l ) - (d„ - d,,) (22) 
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Figure 7. Contact Region with Tangential Loading 


If Pj 
U 


and 

lx 


P 2 are in the non-slip ("stick") region s x is zero, and then 

u * - - d* - 


(23) 


where D x is the global tangential displacement of the cylinders. 

6. NUMERICAL PROCEDURES 


Equations (15) and (18), together with the constraint equations of the foregoing part, 
may be used to determine the modal normal and tangential forces f t and g,. A difficulty 
which arises however, is that the extent of the contact region is not known a-priori. Hence, 
an iterative procedure needs to be developed which will determine the contact region as 
well as the nodal forces. 

For the development of this procedure, it is convenient to superimpose Equations (15) 
and (18) and express them in the form: 


*0x " U ix 


(l-2v)(l+v) 

E 




(24) 
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< 


i j 



and 


.... . (iiMii?) a. ... - tt*k. f . 

“Oz u iz jj u iJ®J £ 1JJ 


(25) 


where c 0 j and d^ are defined as 


e ij - C ij C oj 


and d. - d. • - d - 

1J 1J OJ 


(26) 


where c™ and d™ are influence coefficients of the origin [see Equation (16) and (17).] Then 
by substituting form Equations (24) and (25) into Equation (20) we have 

A.,6.^ + - -bj i,j - 0,...,±n (27) 


where 


*i 


(1-v?) 


+ 


(1-vj) 

*2 


and X 2 


(l-lVjXl+Vj) 


(1-2v 2 )(1+v 2 ) 

e 2 


(28) 


with Ej, j, E 2 , and 2 being the elastic constants of Bodies 1 and 2. 

Similarly, for those nodes in the stick region Equations (24) and (25) together with 
Equations (21) lead to 

Xjcljjf] - X/jjgj - 0 j-0,...,±n (29) 


X 


3 


(1 — 2v 1 )(l+v 1 ) 

E| 


(l-2v 2 )(l + v 2 ) 

E2 


and X 4 


(1-vj) 

E, 


+ 


(1-v^) 

e 2 


(30) 


Finally, for those nodes in the slip region we have 

gi - ±nfj i - -n,...,p,p+m,...,n (31) 

where the sign is chosen so that the force direction is opposite to the direction of slip. 

In Equations (27), (29), and (31) there are 2n+l unknown f s and 2n+l unknown g ( 
(i = 0,...,±n)--that is 4n + 2 equation. Equations (27), (29), and (31) constitute 4n + 2 
equations. However, the system is not complete since when i is zero, Equations (27) and 
(29) are degenerate. Hence, two additional equations are needed. These are obtained from 
global force summation leading to: 
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and G - gjAj 


(32) 


fA 


where Aj is the element area associated with node i. 

The numerical procedure is then: 

(1) Assume a contact region 

(2) Divide the contact region into elements or cells 

(3) Assume a stick regions within the contact region 

(4) Solve Equations (27), (29), and (31) for the nodal forces 

(5) Adjust the extent of the contact region by deleting nodes with negative normal 
forces 

(6) Adjust the stick region using Equation (31) 

(7) Repeat steps (1) to (6) until convergence is obtained 

7. EXAMPLE SOLUTIONS 

To illustrate the efficacy of the method we first considered two frictionless, steel 
cylinders pressed together with a force of 1000 pounds per unit axial length and no 
tangential loading. The cylinders had radii: 100 and 150 inches; elastic modulus: 30 x 10 6 
pounds per square inch; and Poisson ratio: 0.3. The cylinder surfaces were taken to be 
frictionless so that a comparison could be made with the classical Hertz solution, which 
assumes frictionless contact. 

Table 1 shows a comparison of numerical results and the Hertz solution. 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region (in) 

0.136156 

0.137287 

0.83% 

Maximum Normal Pressure (psi) 

9351.35 

9358.90 

0.08% 

Maximum Shear Stress (psi) 

2808.05 

2810.26 

0.08% 

Location of Maximum Shear 
on Z-Axis (in) 

0.05310 

0.05352 

0.79% 


Table 1. Comparison of Hertz and Numerical Solutions for Contacting Steel Cylinders 
Without Friction 

Next to examine the effects of friction (and hence, the consequences of neglecting 
friction), we considered two unlubricated cylinders with radii 100 and 150 inches made of 
steel and aluminum respectively. Different materials were chosen to illustrate the effects 
of material properties on the solutions. (The Hertz solution is based upon identical material 
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properties of the contacting bodies.) As in the first example, the cylinders were pressed 
together with a normal force of 1000 pounds per unit axial length. The elastic moduli for 
steel and aluminum were taken as: 30 x 10 6 and 12 x 10 6 pounds per square inches and the 
Poisson ratios as: 0.3 and 0.33 respectively. The friction coefficient was assigned as: 0.5 
to simulate the unlubricated surfaces. Table 2 shows a comparison of numerical results and 
the Hertz solution. Figure 8 depicts the loading distribution in the contact region. 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region (in) 

0.18896 

0.18873 

0.12% 

Maximum Normal Pressure (psi) 

6738.12 

6810.89 

1.08% 

Maximum Shear Stress (psi) 

2021.43 

2081.32 

(2031.28) 

2.95% 

(0.48%) 

Location of Maximum Shear 
on Z-Axis (in) 

0.07369 

0.07967 

(0.06461) 

8.1% 

(13%) 


Table 2, Comparison of Hertz Solutions (Identical, Frictionless Materials) and 

Numerical Solutions for Contacting Steel/Aluminum Cylinders with Friction. 
(The numbers in parenthesis refer to the aluminum cylinder). 



a = half of contact width from numerical solution 
P 0 = maximum normal pressure from numerical solution 

normal pressure distribution 

□ tangential pressure distribution 

(directions are indicated by arrows) 


Figure 8. Contact Region Loading Distribution for Contacting Steel/ Aluminum 
Cylinders With Friction 
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Finally, to examine the effects of sliding, the steel/aluminum cylinders of the previous 
example were also loaded tangentially to produce sliding. Table 3 shows a comparison of 
the numerical results and the Hertzian solution. 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region (in) 

0.18896 

0.19041 

0.77% 

Maximum Normal Pressure (psi) 

6738.12 

6737.87 

0.00% 

Maximum Shear Stress (psi) 

2021.43 

3365.39 

66.4% 

Location of Maximum Shear 
on Z-Axis (in) 

0.07369 

0.00 

100% 

Shift of Center (in) 

0.00 

0.01123 

100% 


Table 3. Comparison of Hertz Solutions (No Sliding) and Numerical Solutions for 
Contacting/Sliding, Steel/Aluminum Cylinders. 

Figure 9 shows the normal pressure distribution obtained from the Hertz and numerical 
solutions. (Observe the shift in the numerical distribution due to the sliding.) Figures 10 
and 11 show the shear stress distributions on the steel and aluminum cylinders for the Hertz, 
friction and sliding cases. 



a = half of contact width from Hertzian contact 

= maximum normal pressure from Hertzian contact 
□ sliding contact 

Hertzian contact 


Figure 9. Comparison of Hertz and Numerical Solutions for Normal Pressure for 
Contacting/Sliding, Steel/Aluminum Cylinders 

























T maximum shear stress 
P Q maximum norma] pressure (Hertz) 
a half of contact width (Hertz) 


Figure 10. Shear Stress Distributions on the Steel Cylinder. 



T maximum shear stress 
P Q maximum normal pressure (Hertz) 
a half of contact width (Hertz) 

Figure 11. Shear Stress Distributions on the Aluminum Cylinder 
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8. APPLICATION WITH INVOLUTE SPUR GEARS 


To study contact stresses in mating gear teeth we model the teeth as contacting 
cylinders as depicted in Figure 12. The radii of the cylinders depends upon the position of 
the contact point along the line of contact. 

In Figure 12 Q and Q are the centers of the meshing gears and P is the pitch point. 
Ip I 2 , Oj, O z , P, and O are points along the pressure line which is tangent to the base circles 
of the gears as shown. Ij and I 2 are at the intersections of the pressure line and lines 
perpendicular to the pressure line passing through Q and C 2 respectively. O x and 0 2 are 
at the intersections of the pressure line and the addenda of the gears. O is a typical contact 
point along the pressure line. 



Figure 12 Modelling Mating Gear Teeth as Contacting Cylinders 

With this nomenclature , if P t and P 2 are any two points, let PjP 2 be the line segment 
connecting P t and P 2 and let | PjP 2 1 be the length of this line segment. Then the base, pitch, 
and addenda circle radii are 
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(33) 


* 


«b, - IC. 1,1 , r b2 - ICjIjI , r pl - IC.PI , 

r P 2 ” ^ 2 ^ * r «i “ ICjOjI , - IC 2 0 2 l 

where the subscripts 1 and 2 identify the associated gear. 

Since these radii are known for a given pair of meshing gears, it is convenienct to 
express the tooth curvatures (and hence, the contacting cylinder curvatures) in terms of 
these radii. From the properties of the involute profile, the tooth radii of curvature at the 
contact point are simply the distances from the contact point to and I 2 respectively. 
Consider the following contact positions: 

1. Contact at the Pitch Point 

In this case the radii of curvature of the mating tooth surfaces are: 

Pl - Ipl - [Ic^l 2 - lc^l 2 ] 1/2 - ^ - r b 2 ! 

and 

p 2 - Ipl - [lq?l J - IC^I 2 ] 1 ' 2 - - 4 

where the subscripts 1 and 2 identify the associated gear. 

2. Contact on the Addendum Circle of Gear 1 
In this case the radii of curvature are: 

p, - II7o;i - [IC707I J - lc7i;i 2 ] , « - /rf, - r b 2 , 

p 2 - 117071 - 117^1 - H707I - (r pl . r p 2 )sinit> - p, 

where <p is the pressure angle. 

3. Contact on the Addendum Circle of Gear 2 
In this case the radii of curvature are: 

p 2 - \l^0 2 \ - [IC^I 2 - IC^I 2 ] 1 / 2 - - 4 


(34) 

(35) 

(36) 

(37) 

(38) 
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(39) 


* 


and 

Pi - ■Ml 1 - l ^®2 1 “ + " P 2 


4. Contact at Other Points 

Cases 2 and 3 are at the contact extremes with Case 1 being intermediate. For some 
other contact point, the radii of curvature will be different than these cases, but it can be 
calculated using the same procedures as in the above cases. 


The developed numerical procedures may now be used to calculate the contact 
stresses on the mating teeth. To illustrate this applicaiton, consider two identical spur gears 
in mesh with the following properties. 

Number of teeth: 18 

Pitch diameter: 3.5433 in. 

Face width: 0.3937 in. 


Pressure angle: 

Transmitted load: 
(tangential) 

Elastic modulus: 

Poisson ratio: 


20 deg. 

808 lb. 

30 x 10 6 psi 
0.3 


This data yields the following parameters: 

Base circle radius: 1.6648 in. 

Addendum circle radius: 1.9685 in. 

Contact force: 860 lb. 

(normal) 

Contact force per 

unit length: 2184 lb. 

To put the results in perspective, it is informative to make comparisons between the 
maximum contact stresses and the maximum fillet (root) stresses, which are documented 
extensively in the technical literature. 
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Table 4. presents the results and comparisons. The tensile fillet stresses are obtained 
using the Lewis formula [48], the Lewis factor [49], and the stress concentration factor [50, 
51]. The compressive fillet stress may be obtained using finite element methods [52]. 
Finally, the stresses of Table 4 are listed in psi. The maximum Von Mises stresses occur 
beneath the surface at depths of 0.0056 in for pitch point contact and 0.0038 in for 
addendum point contact. The maximum compressive contact stresses occur on the surface. 



Pitch Point Contact 

Addendum Contact 

Fillet Tensile Stress 

N/A 

46,217 

Fillet Compressive Stress 

N/A 

55,080 

Compressive Contact Stress 

194,480 

286,013 

Von Mises Contact Stress 

58,407 

85,899 


Table 4. Stress Results and Comparisons 


9. DISCUSSION 

Although these examples employ bodies with simple geometries (cylinders), the method 
outlined herein is not restricted to simple geometries. Indeed, there are no restrictions on 
the geometry so long as the contacting surfaces are continuous and geometrically smooth 
(that is, with continuous derivatives). 

The examples demonstrate the efficacy of the method. They show that the numerical 
results are extremely close to those of the classical Hertz solutions, for those cases where 
the Hertz solution is applicable. Indeed, the stresses are virtually identical whereas there 
is a slight difference in the geometric results - that is, the contact region width and the 
location of the maximum shear stress. (The Hertz solution assumes frictionless static 
contact between identical materials.) 

The second and third example demonstrate the significant effects of friction and sliding 
in contact analysis. The examples also show that the softer material has lesser stress. 

In addition it is seen that friction and sliding significantly increase the shear stresses. 
Indeed, when there is sliding with a relatively high friction coefficient, the maximum shear 
stress occurs on the surface as opposed to being beneath the surface (See Figures 10 and 
11 ). 

Finally, in the examples with the gear stresses it is seen that the fillet stress is dependent 
upon the loading and geometry whereas the contact stress depends upon the elastic 
properties as well. 

Regarding the method itself it is seen that it employs a double discretation: one of the 
contact region, the other of the force distribution. The execution of the method involves 
an iterative procedure to determine the extent of the contact region. In this way the method 


♦ 
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is clearly more laborious than the classical Hertz method and even elementary finite- 
element methods (FEM). In exchange, however, the method is believed to provide a more 
accurate and comprehensive analysis than either the Hertz or FEM techniques. 

The method outlined herein is developed for a two-dimensional analysis with a 
triangular load on the individual elements or cells. The triangular load provides for a 
piecewise linear representation of the contact loading. The extension to higher order load 
representations (for example, the use of cubic splines), and to three-dimensional analysis, 
is conceptually "straight-forward." Such extensions are currently being developed. 

10. CONCLUSIONS 

A point load superposition procedure based on the theory of elasticity has been used 
to study contact stresses between elastic bodies, including meshing spur gear teeth. This 
method is validated by comparison with results from the classical Hertz method for 
frictionless cylinders in contact. The following specific conclusions are reached: 

1. A new method of contact stress analysis has been presented. It is a numerical method 
which has fewer restrictions than classical procedures. 

2. The method is applicable for a broad range of geometries of the contacting bodies. The 
contacting bodies may have different material properties. The contacting surfaces need not 
be smooth. Sliding between the contacting surfaces is also permitted. 

3. The accuracy of the method is demonstrated by a comparison with results from classical 
methods for simple cases where the classical method is applicable. 

4. The established efficacy of the method justifies applications with contacting gear teeth. 
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